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1.  Summary  of  the  Accuracy  Test.  The  valuea  in  the  parentheses  are  the 

modified  percent  errors  where  a  single  extremum  percent  error  at  the  zero 

crossing  is  eliminated  from  each  case.  Other  values  are  Si  =  2.0, 

a§  =  0.1,  yj  =  -0. 1,  y2  =  0. 1,  b'  =  1/412,  and  «  for  the  termination 

criterion  is  10" 6  32 


Unified  Real  Part  of  the  Susceptibility  for 

Collisional  Broadening 


1.  INTRODUCTION 


As  a  consequence  of  the  extensive  study  of  the  emission  and  absorption  of  spec- 

tral  radiation,  significant  advances  have  been  made  in  the  understanding  of  the  imagin- 
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ary  part  of  the  susceptibility  for  collisionally  broadened  transitions.  An  absorption 
line  shape  in  the  impact  approximation  has  been  developed  that  is  applicable  from  the 
microwave  through  the  infrared  and  for  which  a  first  order  treatment  of  line  coupling 
is  available.  In  this  report  we  develop  the  real  part  of  the  susceptibility,  applicable  to 
the  computation  of  the  index  of  refraction  and  propagation  delays,  which  is  associated 
with  this  absorption  line  shape. 

(Received  for  publication  12  August  1986) 

1.  Clough,  S.  A.,  Davies,  R.W.,  and  Tipping,  R.H.  (1983)  The  line  shape  for  collision - 

ally  broadened  molecular  transitions:  a  quantum  theory  satisfying  the  fluctuation 
dissipation  theorem,  in  Spectral  Line  Shape  Vol.  2,  edited  by  K.  Burnett, 

Walter  de  Gruyter,  New  Vork,  NY,  1?83,  pp  553-568. 

2.  Liebe,  H.J.  (1981)  Modeling  attenuation  and  phase  of  radio  waves  in  air  at  frequen¬ 

cies  below  1000  GHz,  Radio  Science,  16,  No.  6:1183-1199. 

3.  Kemp,  A.J.,  Birch,  J.R.,  and  Afsar,  M.N.  (1978)  The  refractive  index  of  water 

vapor:  a  comparison  of  measurement  and  theory.  Infrared  Physics,  18:827-833. 

— —  '■  '  "  VW 

4.  Marshall,  A.G.,  and  Roe,  D.C.  (1978)  Dispersion  versus  absorption:  spectral  line 

shape  analysis  for  radiofrequency  and  microwave  spectrometry,  Analyt.  Chem., 
50,  No.  6:756-763.  - 


The  complex  susceptibility,  X  ,  characterizes  the  linear  relationship  between  the 
polarization,  P,  and  the  impressed  field,  E  =  EQexp(i2irc >>t) ,  through  the  expression 

P  =  X  ( v )  EQexp(i2*c  vt) ,  (1.1 

where  v  is  the  wavenumber  value  and  c  is  the  velocity  of  light.  The  susceptibility  is 
expressed  in  terms  of  its  real  and  imaginary  parts  as 

X  (i»)  =  X'(v)  -  i  X"(v).  (1.2 

For  a  plane  wave  propagating  in  the  z  direction  through  a  nonmagnetic  dielectric  medi 
urn  (permeability  =  1)  with  complex  index  of  refraction  n , 


we  have -for  the  electric  field, 

E^  =  E0exp[-2»  v  n"z]  exp[i2ir  v  (ct-n'z)] ,  (1.4 

for  the  magnetic  field, 

Hy  =  n  Egexp[-2 *  v  n"z]  exp[i2wK  (ct-n'z)] ,  (1.5 

and  for  the  real  part  of  the  cycle-averaged  Poynting  vector, <N>,  describing  the  spa 
tial  distribution  of  the  energy  density, 


1  0  2 

<  N  >  =  —  ( - )  n'  E«  exp[-4*  v  n"] . 

z  2  4»  0 


As  established  in  Appendix  A,  this  leads  to  relationships  between  the  susceptibility 
and  the  index  of  refraction 


=  n’2  -  1  -  n"2, 


4»X"  =  2  n'  n", 


and  to  the  definition  of  the  absorption  coefficient  per  unit  length  (1/cm),  *  , 


For  gases,  Eqs.  (1.7),  and  (1.8)  are  often  expressed  as 


4  *  X'  =  2(n'  -1) ,  (1.7a) 

and 

4  #  X  ”  =  2n" ,  (1.8a) 

2  2 

where  n'  has  been  taken  as  =:  1,  and  (n*  -  1)  >>  n"  .  Therefore,  the  absorption  co¬ 
efficient  for  gases  is  given  in  terms  of  the  susceptibility  by 

<c  =  8  ir2  v  X" .  (1.10) 

The  expression  for  the  imaginary  part  of  the  susceptibility  in  terms  of  the  power 
spectral  density  function  <t>  ( v )  may  be  developed  in  three  forms, 

(1-e -*") 

X"(«0  =  -  [♦(*)  +  ♦(-»»)],  (1.11a) 

(1V^) 

X"(V  )  =  (l-e'^")  <t>  (  r),  (1.11b) 

and 

X"  (V)  =  <t>(v)  -  <t>  (-v),  (1.11c) 

where  0  =  hc/kT.  These  three  forms  for  X”(v)  are  equivalent  if  <t>  (  v )  satisfies  an 

5  6 

important  theorem,  the  fluctuation  dissipation  theorem  (detailed  balance)  ’  which  is 
given  by 

4>(-v)  =  e  ^v0(i').  (1.12) 


5.  Kubo,  R.  The  fluctuation-dissipation  theorem,  in  Reports  on  Progress  in  Physics , 

Vol.  XXIX  Part  1,  (1966),  A.C.  Stickland,  Ed. ,  Tne  Institute  of  Physics  and  the 
Physical  Society,  London,  pp  255-284. 

6.  Callen,  H.B.,  and  Welton,  T.A.  (1951)  Irreversibility  and  generalized  noise,  Phys. 

Rev.  83:34-40. 


If  ♦  O' )  does  not  satisfy  detailed  balance,  then  the  form  given  by  Eq.  (1.11a)  has  cer¬ 
tain  attractive  properties:  radiation  balance  between  emission  and  absorption  is  main¬ 
tained  independent  of  the  correctness  of  $  ( * ) ,  and  [$(»>)  +  ♦(-•')]  is  an  even  function 
of  v  and  satisfies  important  summation  properties.  ’  It  will  prove  convenient  to  write 
Eq.  (1.11a)  in  the  form 

X"(iO  =  tanh(£*/2)  [*(•/)  +  4(-v)],  (1.13) 

We  now  assume  that  the  imaginary  part  of  the  susceptibility  may  be  obtained  from 
the  superposition  of  spectral  contributions  from  individual  transitions, 

X"(*0  =  £x".(*)  =tanh(/W2)  !.[*.(  v)  +  *.(-  v)] ,  (1.14) 

where  $.(-v)  =  e  ^.(  v ) . 

At  this  stage,  we  make  the  impact  approximation,  obtaining 

X"(*/)  =  tanh(0W2)  2  S.(mp,T)[f.(  v  ,  mp,T)  +  f.(  -  *  ,mp,T)] ,  (1.15) 

3  3 

where  for  radiator  density  m  (molec/cm  ),  perturber  density  m  (molec/cm  ), 

~  1  "  P 

S.(m  ,T)  (cm  )  is  the  intensity  associated  with  the  i-th  transition  and  f.(  v  ,m  ,T)  is 

1  r  *  o  P  9 

the  Lorentz  function.  The  intensity  Sj  in  terms  of  the  transition  strength  p.  (debye  ) 
is  given  by 

S  (m  ,T)  =  m  ,JL_  10*36  P. 2  [p."(T)  +  P.'(T)],  (1.16) 

1  r  3hc  11  1 

for  which  P and  p/  are  the  fractional  populations  of  the  lower  and  upper  states  for 
the  transition,  respectively.  The  Lorentz  function  is  defined  through  the  expression 

1  a.(mp,T) 

f.(,',m  ,T)  =  (— )  -  ,  (1.17) 

1  P  "  2  9 

(•'-•'j)  +  a.^mp.T) 

in  which  v  is  the  transition  wavenumber  value  and  Q.(m  ,T)  is  the  halfwidth.  In  col- 
1  i  P 

lisional  broadening  of  gases,  the  halfwidth  is  a  linear  function  of  perturber  density  ex¬ 
pressed  through  the  relationship 

7.  Van  Vleck,  J.H.,  and  Huber,  D.L.  (1977)  Absorption,  emission,  and  linebreadths: 
a  semihistorical  perspective,  Rev.  Mod.  Phys.  49;  No.  4:939-959. 


(1.18) 


«j(m  ,T)  =  ( - )aio(T) 

mo 


The  integral  property  of  the  function  f.(  v  •  n>p,T), 


y  f.(  v  ,rn  ,T)  d  v  =  1, 

—  oo  1  P 


(1.19) 


preserves  the  summation  property  of  [♦(•>)  +  ^(-«>)].  In  the  impact  approximation , 
the  equivalency  of  Eqs.  (1.11a)  through  (1.11c)  is  not  preserved  since 


f(-  v)  *  e^U  v). 


(1.20) 


For  clarity  of  presentation,  the  explicit  dependence  of  S.(mr,T),  f.(|'.,m  ,T),  and 
o  j(n)p,T)  on  mr.nip  and  T  is  suppressed,  giving  for  the  imaginary  part  of  the  suscept¬ 
ibility 


X"(i»)  =  tanh(/?p/2)  IS.  (f.(»)  +  f.(-i*)]. 
i 


(1.21) 


It  may  readily  be  established  that  Eq.  (1.21)  reduces  to  the  Van  Vleck-Weisskopf  re¬ 
sult  for  flv  <scl  and  to  the  infrared  result  for  ffv  »l. 

Line  coupling  has  become  of  increased  interest  because  of  its  effects  on  spectral 


properties  directly  related  to  the  remote  sensing  problem;  the  60  GHz  band  for  oxygen 

- J  ,1 _  err  r - -  "  1  \  A  1 - V.  f - - U - ri  ^  |n  ...a  ..till  mnlnHo  tha 


and  the  o67  (cm  )  Q  branch  for  carbon  dioxide.  In  this  report,  we  will  include  the 
effects  of  a  first  order  treatment  of  line  coupling  that  also  satisfies  the  condition  of 
superposition.  We  have  for  the  line  coupling  case  in  the  impact  approximation 


X"(  v )  =  tanh(/?v/2)  I  S.(fj(  v  )  +  y.g.(  v  )  +  f.(-  v)  +  y.g.(-  O] , 


(1.22) 


where  the  line  coupling  coefficient  y^  is  linearly  dependent  on  pressure  and  the  func¬ 
tion  g.(  v  )  is  defined  as 


gjOOM-) 


(v-  v.) 


(1.23) 


(v  -  v  ,)2  +  «  } 

l  l 


8.  Van  Vleck,  J.H.,  and  Weisskopf,  V.F.  (1945)  On  the  shape  of  collision-broadened 
lines.  Rev.  Mod.  Phys.  17,  Nos.  2  and  3:227-236. 
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We  will  consider  the  real  part  of  the  susceptibility  in  the  impact  approximator 
given  by  Eq.  (1.21),  for  two  limiting  spectra)  domains,  the  millimeter  and  infrared 
regions,  for  which  analytical  results  are  available.  A  proposed  analytical  approximating 
function,  consistent  with  the  results  in  these  two  limiting  spectral  regions,  is  then  devei 
oped  for  the  real  part  of  the  susceptibility.  We  also  develop  an  expression  for  the  real 
part  of  the  susceptibility  by  taking  the  Hilbert  transform  of  Eq.  (1.21)  in  which  the  tanh 
function  is  represented  as  a  series  expansion  in  terms  of  Lorentz  functions.  This  result 
is  then  evaluated  numerically  to  establish  the  validity  of  the  proposed  approximate  func 
tion  over  the  entire  spectral  domain. 

Since  the  impact  approximation  does  not  provide  an  appropriate  description  for  the 
line  wing,  as  evidenced  by  its  failure  to  satisfy  detailed  balance  and  to  decay  exponen 
tially  in  the  wing,  we  will  discuss  in  Section  5  the  effect  of  the  wing  contribution  to  the 
real  part  of  the  susceptibility. 

2.  THE  KRAMERS-KRONIG  RELATIONSHIP  AND  THE  HILBERT  TRANSFORM 

q_  i  g 

The  Kramers-Kronig  relationship  plays  the  central  role  in  our  development  by 
relating  the  real  part,  X'(v),  and  the  imaginary  part,  X”(  v),  of  the  complex  suscept¬ 
ibility,  X(v),  through  the  Hilbert  transform  as 

r°°  1 

X'(»)  =  H[  X  "(  v  )]  =  ( 1/  n  )  /[ - ]  X  "(u)  du,  (2.1a) 

--  (u  -  v  ) 


9.  Bolton,  H.C.  (1969)  Some  practical  properties  of  Kronig-Kramers  transforms.  The 
Philosophical  Magazine,  19,  No.  159:487-499. 

10.  Fahrenfort,  J.  (1963)  The  methods  and  results  of  dispersion  studies.  Chapter  XI 

in  Infra-red  Spectroscopy  and  Molecular  Structure,  M.M.  Davies  Ed.,  Elsevier, 
New  York,  NY. 

11.  Huber,  D.L.,  and  Van  Vleck,  J.H.  (1966)  The  role  of  Boltzmann  factors  in  line 

shape,  Rev.  Mod.  Phys.  38,  No.  1:187-204. 

12.  Kramers,  H.A.  (1929)  Die  dispersion  und  absorption  von  Roentgenstrahlen , 

Physikalishe  Zeitschrift,  30:522-523. 

-  VW 

13.  Kronig,  R.  de  L,  and  Kramers,  H.A.  (  1928)  Absorption  and  dispersion  in  X-ray 

spectra,  Z.  Physik,  48,  No.  3-4:174-179. 

14.  Kronig,  R.  de  L.  (1926)  On  the  theory  of  dispersion  of  X-rays,  J.  Opt.  Soc.  Am. 

12,  No.  6:547-557.  - , 

15.  Martin,  P.C.  (1967)  Sum  rules,  Kramers-Kronig  relations,  and  transport  coeffi¬ 

cients  in  charged  systems,  Phys.  Rev. ,  16J,  No.  1:161-155. 

16.  Toll,  J.S.  (1956)  Causality  and  the  dispersion  relation:  logical  foundations, 

Phys.  Rev.  104,  No.  6:1760-1770. 


and  the  inverse  transform  as 


Xn(.v )  =  H'^X  '(v)]  =-(!/«) 


r°°  1 

J  [ - j  x  .(u)  du. 

•a  (U  -  V  ) 


(2.1b) 


It  is  noted  that  the  Hilbert  transform  is  unique  up  to  an  additive  constant  and,  there¬ 
fore,  is  defined  for  functions  that  vanish  at  infinity. 

17 

This  relationship  has  been  expressed  in  several  alternative  forms  including 
those  described  below: 


X'C*' )  =  (2 /»  ) 


r°°  1 

/  I - ]  u  X"  (u)  du, 

o  (u2-,2) 


X  "(V)  =  -  (2/ir  ) 


r 


(u2- 


-]  v  X’(u)  du, 


w2) 


(2.2a,  Ref.  10) 


(2.2b,  Ref.  5) 


or 


X’(  v)  =  (2  /  »  ) 


■]  (u  X"(u)  -  »>  X  "(v)]  du, 


X"(x)  =  -  (2 In  ) 


■] 


[  X  '(u)  -  X  '(  v )]  du. 


(2.3a,  Ref.  15) 


(2.3b,  Ref.  17) 


The  equivalence  of  Eqs.  (2.1)  and  (2.2)  can  be  shown  by  breaking  the  integral  in  Eq. 
(2.1)  into  two  parts  for  the  negative  and  positive  range  of  the  argument  u  and  recom¬ 
bining  them  using  the  oddness  of  x"(  v  ) ,  or  the  evenness  of  X '(  v ) .  On  the  other 
hand,  Eqs.  (2.2)  and  (2.3)  are  equivalent  since  the  principal  value  of 


P 


/ 


Q 


(ii 


*2) 


du  =  0, 


(2.4) 


for  any  constant  Q. 


17.  Townes,  C.H.,  and  Schwalow,  A.L.  (1955)  Microwave  Spectroscopy,  McGraw-Hill 
Book  Company,  New  York,  NY. 


The  relationship  of  the  real  and  imaginary  parts  of  the  susceptibility  can  also  be 
discussed  in  terms  of  the  Fourier  transform.  Let  Y(t)  be  the  inverse  Fourier  transform 
of  the  complex  susceptibility  X(f),  i.e. ,  Y(t)  =  7  1  [  X  ( f ) ] .  Note  that  for  consistency 
with  the  conventional  definition  of  the  Fourier  transform ,  x  ( . )  is  a  function  of  f=  »  c 
with  the  frequency  f  in  Hertz,  whereas  the  Fourier  transform  pair  is  defined  in  terms 
of  time  t  (sec)  and  frequency  f  (Hz)  as  below. 


G(f)  =  7  [g(t)l  =  /g(t)e"i2  '  ft  dt, 


(2.5a) 


g(t)  =  7'1[G(f)]  =  y*G(f)e+l2  *  ft  df.  (2.5b) 

-  OB 

Now,  Y(t)  is  clearly  a  real  and  causal  function  because  it  represents  a  causal  process 
of  the  wave-medium  interaction.  Therefore,  its  Fourier  cosine  and  sine  transforms  can 
be  associated  with  X,(»)  and  X"(v )  respectively  as  follows: 

oo 

X  (f)  =  7  [Y(t)l  =  f  e  12  ”  ft  Y(t)  dt. 


=  f  (cos2»  ft  -  i  sin 2  it  ft)  Y(t)  dt, 

O 

x(f)  =  (1/2)  7clY(t)J  -  i  (1/2)  7s[Y(t)l, 

=  X  '(  v)  -  i  X"(  v),  (2.6) 


where  ^7  c [ .  ]  and  7  a t .  ]  denote  the  Fourier  cosine  and  sine  transforms  defined  as 


Gc(f)  =  7c[g(t)]  =  2  y"cos2  irft  g(t)  dt, 

O 

g(t)  =  7  c“  1[G(f)J  =  2  /« js2  eft  Gc(f)  df. 


(2.7) 


and 


Gs(f)  =  7  s[g(t)l  =  2 
g(t)  =  7g'1[G(f)]  =  2 


(2.8) 


Since  for  the  causal  function  Y(t),  Y(t)  =  0  for  t  <0. 


Hence  we  have 


X'(f)  =  (1/2)  7cIY(t)]. 


x"(f)  =  (1/2)  7  (Y(t)] . 

O 


Note  that  the  change  of  variable  is  also  needed  here.  After  some  algebra,  the  following 
alternate  relationship  can  be  obtained. 


X'(f)  =(1/4)  7  [  7  " 7  X  "( v )]] , 


(2.10) 


or,  in  an  equivalent  form 


C(f)  =  /  cos2  *  ft  f  si 

0 


sin2»  ut  X"(u)  du 


(2.11) 


It  is  noted  that  all  the  real  and  imaginary  parts  appearing  in  these  expressions 
have  the  required  symmetry  properties;  the  real  parts  are  even  functions  and  the  im¬ 
aginary  parts  are  odd  functions. 

S.  IMPACT  LINE  SHAPE 

We  consider  the  single  transition  susceptibility  in  this  section  for  the  sake  of  no- 
tational  simplicity.  In  Eqs.  (1.21)  and  (1.17),  we  eliminate  the  summation,  utilize  v 
in  place  of  Vj,  and  drop  all  other  i  subscripts.  These  may  be  reinvoked  as  the  need 
arises. 

The  imaginary  part  of  the  complex  susceptibility  at  wavenumber  value  v  due  to  a 
transition  at  with  intensity  S  (cm  S  and  half- width  «(cm  *)  can  be  expressed  as11 

X"(x)  =  tanh(/3 v  12)  S  [f(v)  +  f(-  iO],  (3.1) 

where  f(  v)  is  the  Lorentz  line  shape  function  as  before. 


f(„)=-L - 1 


(1/  -  V  0)2  +  <*  2 
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It  is  useful  to  define  a  related  function  g(  v)  where 


£(*')  =  — 


0> 


c'  -  V  + 


(3.3) 


and  to  establish  the  following  identities  involving  f(  »)  and  g(  v)  that  will  be  conveni¬ 
ent  for  this  study. 


First  we  consider  the  millimeter  wave  region  for  which  the  impact  result  for  X  ( v ) 


8 


is  well  known.  For  this  spectral  region,  with  fi  v  «■  1,  we  take 
tanh(£  i//2)<»  0  v  12, 
giving  for  the  imaginary  part. 


(3.5) 


V"(  v )  :mm  =  (/8  »/2)S  [f(v)  +  f(-v)]. 


(3.6) 


Then  the  corresponding  real  part  x’(  vi  :mm  is  given  by 


X'<  v):ron)  =  HI  x"(  *’)  -mm]  =  (S0/2)HI»*f(v)  +  vf(-  „)) 


=  (S  [v0f(  r)  +  «g(  v)]  -  [  vQf(- v)  +  «g(- v)  J  . 


(3.7) 


The  Hilbert  transforms  appearing  in  this  equation  are  slight  modifications  of  stand- 
18 

ard  ones,  but  are  derived  in  the  appendix  because  of  their  significance  to  our  develop¬ 


ment.  Using  these  results,  we  arrive  at  the  final  result  for  the  millimeter  region. 


18.  Bateman  Manuscript  Project,  (1954)  Tables  of  Integral  Transforms,  Vols.  1  and  2, 
McGraw-Hill,  New  York,  NY.  - 
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vf(v)  =  vQf(  v)  +  «g(  v). 

(3.4a) 

Vg(  K)=-«f(v)+  VQg(  v)  +  1  /»  , 

(3.4b) 

i 

vf(  k  )  =  -  *'„«-»')  -  «g(-v), 

(3.4c) 

•'g(-v)  =  af(-v)  -  v0g(-v)  -  1/s. 

(3.  4d) 

$ 

These  can  be  shown  straightforwardly. 

$ 

>> 

>1 


>;i 


4 


i 


I 

3 


.V 


« 


In  the  infrared  case,  for  which  (3 v  12  »1,  tanh(0  w  12)  may  be  approximated  by 
unity.  Since  the  Hilbert  transform  is  defined  over  the  entire  real  line  and  since  it  is 
essential  to  keep  the  symmetry  property  of  tanh( . ) ,  we  utilize  the  approximation 

tanh(0v/2)=  sgn(v),  (3.9) 

where  the  signum  function  sgn(  v)  takes  the  value  1( - 1)  if  the  argument  v  is  positive 
(negative).  Then  X"(v)  can  be  approximated  by 

Y"(v):m  =S  sgn(„)  [f(v)  +  f(- v)]  .  (3.10) 

Then  the  corresponding  real  part  of  the  susceptibility  is 


V'(v):m  =H[Y  ”(v):IR]  =  (l/<r  )  f  x"(u) 


IR, 


(u- v)  du. 


(3.11) 


Evaluating  the  principal  value  of  the  integrals  in  this  equation ,  we  obtain , 

-1, 


X'(v):m  =  -  (2S/  *  )  tan'^Vp/e)  [g( *)  +  g(-v)] 

2  2 
%  +  « 

+  (S /»  )  In  [ - 1  [f(  v)  +  f(-  *)]  . 


(3.12) 


Using  the  infrared  approximation  Vp/a  »1,  and  assuming  v  =  Vp,  this  equation  may 
be  simplified  to 


x’(v):m  =  -  S  [g(*)  +  g(- v)]  .  (3.13) 

The  elimination  of  the  logarithm  and  the  setting  of  the  tan  1  function  at  » /2  may  be 
viewed  as  eliminating  the  spurious  discontinuity  of  Eq.  (3. 10)  at  the  origin  (  v  =  0). 

The  discontinuity  is  spurious  since  it  arises  from  approximating  tanh(0  » 12)  by  sgn(  v) 
and  this  approximation  is  not  valid  in  the  vicinity  of  the  origin. 

The  result  given  by  Eq.  (3.13)  for  X'(*’):IRmay  be  obtained  directly  by  approxi¬ 
mating  x"(x)  in  Eq.  (3.11), 


V"(v):IR2  =  s  [fOO  -  t (-•')]  ■ 


(3.14) 
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It  is  illuminating  to  note  that  this  form  of  X  "(  v)  is  consistent  with  the  form  given  by 
one  of  the  alternate  expressions  for  the  imaginary  part  of  the  susceptibility  in  terms  of 
the  power  spectral  density  function,  Eq.  (1.4c).  This  approximate  form  for  X"(v)  has 
an  advantage  over  that  given  in  Eq.  (3.10)  since  it  has  no  discontinuity  at  the  origin, 
and  consequently  will  prove  to  be  useful  when  we  consider  the  line  coupling  problem  in 
Section  4. 

To  study  the  transition  of  X  '(»)  from  one  spectral  regime  to  the  other,  we  first 
consider  the  delta  line  shape  function.  With  the  delta  line  shape  function,  the  exact 
Hilbert  transform  can  be  obtained  without  using  approximations  to  the  tanh(  i  v  12)  func¬ 
tion.  Thus  by  replacing  f(*  )  in  Eq.  (3. 1)  by  the  delta  function  «  (  v  -  v  0) ,  we  obtain 

X  '(•')=  H  [  tanh{0  *  /2)  S  [  4  (  *  -  **0)  +  4  (  v  +  kq)]] 


1  1 

=  S  tanh(/J««  .12)  [  -  +  - 1 

(-O'  *>  («'0+  •'> 


=  2S  tanh(0 *>q/2) 


(3.15) 


This  expression  suggests  that  if  the  line  shape  function,  f(  r).  is  narrow  compared  to 
the  variation  of  the  tanh(0  v  12)  function,  then  the  transition  of  the  magnitude  of 
x  '(*')  from  the  millimeter  to  the  infrared  is  specified  by  the  same  tanh(  0  v  12)  func¬ 
tion  as  for  X  "(f)  except  that  tanh(0  v  12)  is  evaluated  at  the  line  center  v  =  v  q. 

Considering  the  two  results  in  Eqs.  (3.8  and  (3.13),  and  the  above  observation, 
we  now  postulate  the  following  real  part  of  the  susceptibility  for  use  in  both  millimeter 
and  infrared  regions. 

X  ’  (v)  =  S  tanh(0i»  fl/2)  [(  gd-)  +  (  »/>'„)«»')]] 

+  [-g(-x)  +  (<*/i»0)f(-i>)]]  .  (3.16) 

This  result  reduces  to  the  millimeter  result  in  Eq.  (3.8)  by  letting  &  *'q/2  become  small, 

and  to  the  infrared  result  in  Eq.  (3.13)  by  letting  0  vQ/2  become  large  and  v  small. 

We  now  perform  the  inverse  Hilbert  transform  of  X  '  ( i> )  to  obtain  X  "  (  v  )  for 

P  P 

comparisons  with  previously  used  expressions, 

X  "p(0  =  H1[X'p(»')]  =  -H  [X'p(i')] 

=  S  (  v/vQ)  tanh(  0i»q/2)  [  f(  v  )  +  f(-i/)  ]  . 


(3.17) 


i«iiii*iMiimmnuiwiwun»uPuwwuuu»uuwn 


For  the  millimeter  case,  the  approximation  tanh(  P  v  12)  a.  Pv  12  reduces  this  result  to 
Eq.  (3.6).  On  the  other  hand,  for  the  infrared  case,  tanh(  p  1^/2)  —  1.  Eq.  (3.17) 
gives 

X"  (•'):m  =  S  (  v!  i/Q)  [f(y)  +  f(- v)]  .  (3.18) 

The  expression  states  that  the  proposed  real  part  gives  rise  to  a  linear  approximation 
v  I v  q  f or  the  tanh(  P  v/2)  function  in  the  infrared  region.  It  is  essentia]  to  note  that 
this  function  is  a  valid  approximation  in  the  infrared  for  frequencies  near  the  line  cen¬ 
ter  for  which  v  is  large  and  v— 

The  tanh(  p  *'q/2)  dependence  in  Eq.  (3.16)  suggests  that  the  proposed  real  part 
might  be  obtained  through  the  line  center  approximation  defined  by  the  replacement  of 
tanh(  P  v  12)  by  tanh(  p  v^/2)  and  then  taking  the  Hilbert  transform.  However,  the 
existence  of  the  second  and  the  fourth  terms  in  x  ’  (  v)  proves  otherwise.  It  will  be 
shown  later,  however,  that  the  major  contribution  to  the  Hilbert  transform  comes  from 
the  line  center  rather  than  the  line  wing.  This  gives  a  rationalization  for  the 
tanh(  P  v 0/2)  dependence  of  the  proposed  real  part. 

4.  IMPACT  LINE  SHAPE  WITH  LINE  COUPLING 


Although  the  inclusion  of  line  coupling  in  the  determination  of  the  imaginary  part 

of  the  susceptibility  for  a  molecular  system  is  in  general  a  complex  problem ,  a  first 

19-21 

order  treatment  has  been  developed  in  which  superposition  is  preserved.  In  our 

formalism ,  the  imaginary  part  of  the  susceptibility  may  be  expressed  as , 


■'VI 


,(x)  =  I  Sj  tanh(  P  » 12)  £f.(  v )  +  yjg.(  v  ) 


+  v- 


v  )  + 


Vigjf-O] 


(4.1) 


where  yj  is  the  line  coupling  coefficient,  which  in  general  will  be  subject  to  the  con 
straint2* , 


1 

i 


Vi 


=  0  , 


(4.2) 


19.  Baranger,  M.  (1958)  Problem  of  overlapping  lines  in  the  theory  of  pressure 

broadening,  Phys.  Rev.  Ill,  No.  2:494-504. 

20.  Rosenkranz,  P.W.  (1975)  Shape  of  the  5  mm  oxygen  band  in  the  atmosphere, 

IEEE  Trans.  Ant.  Propagat.,  AP-23,  No.  4:498-505. 

21.  Smith,  E.W.  (1981)  Absorption  and  dispersion  in  the  Oo  microwave  spectrum  at 

atmospheric  pressures,  J.  Chem.  Phys. ,  74:6658-6673. 
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referred  to  as  the  zero- sum  rule.  The  spectral  properties  of  X  "( »<)  for  overlapped 
lines  are  dependent  on  this  zero-sum  rule,  and  the  relationship  of  the  y  .  to  the  a  . . 

At  this  stage  in  the  development,  we  drop  the  summation  from  Eq.  (4.1)  and 
adopting  the  notation  of  Section  3,  we  obtain  the  imaginary  part  of  the  susceptibility 
associated  with  a  single  transition, 

X  "c(  v  )  =  S  tanh(  0V  I 2)  (  f(  v  )  +  yg(  v  )  +  f(-  v)  +  yg(  -  v  )  ]  ,  (4.3) 

We  again  start  in  the  millimeter  spectral  regime,  where  tanh(0  v  12)  =>  0  v  12,  for 
which  we  obtain, 


X"c(0:mnl  =  S  (/J„/2)  1  f(a)  +  yg(  v)  +  f  ( -  r  )  +  ygt-r)  ]  .  (4.4) 

Using  the  previously  established  identities  among  f(  v)  and  g(  v)  given  in  Eq.  (3.4), 
we  obtain 

x"0(  "):mm  =  (S/J/2)  [((i»0  +  ya)f(v)  +  (a  ♦  y  t-0)g(  v  )  ]  +  y/w 

-  |  (»0  -  ya)f(-»)  +  (a  +yc())g(-r)  )  -  y/i  (4.5) 

We  may  cancel  ±y/w  terms  in  this  expression  and  obtain, 

X  V^'inm  =  (S/»/2)  [  t(  vQ  -  ya)f(v)  +  (a  +  y  »'0>g(>'>) 

-  (  ( i>0  -  ya)f(- v)  +  (  a  +  y  »* Q)g< -  v)  1  ]  .  (4.6) 

We  now  take  the  Hilbert  transform  of  Eq.  (4.6)  to  obtain  the  real  part  of  the  sus¬ 
ceptibility  for  the  line  coupled  case  in  the  millimeter  region,  obtaining 


=  (S/J/2)  [l  -  (»0  -  y«)g(»)  +(a+yi>0)f(«')] 

+  l  '  ( *»  q  -  ya)g(-t»)  +  (a  +  y«'jj)f(-i')l  ].  (4.7a) 

X  c(  v  ):mm  =  <S^0/2)  C  1  ‘  (1  '  ya/*'0>K(,'>  +  (a  /  v  0  +  y)f< »» )  1 

+  l-(l-ya/i»0)g(-i')+(a/»'0  +  y)f(-i')l].  (4.7b) 
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Thus,  the  significance  of  Eq.  (4.6)  for  X  "c(  v )  is  that  its  Hilbert  transform  is  well 
defined  and  the  inverse  Hilbert  transform  of  the  corresponding  X  ’  (  v )  matches  the 
original  expression.  Thus  the  difficulty  encountered  in  taking  the  Hilbert  transform  of 
Eq.  (4.4)  in  earlier  studies  is  eliminated. 

It  is  straightforward  to  check  the  agreement  of  the  inverse  Hilbert  transform  of 
this  expression  with  X  "  (  v )  given  in  Eq.  (4.6). 

For  the  infrared  region,  in  which  tanh(£»>/2)  =■  1,  we  use  the  form  given  by 
Eq.  (3.14)  for  the  line-coupled  imaginary  part, 


x’y’O'lR  =S  [|f(*>)  +yg(*)l  -  (H-a)  +  yg(- «')]  ]  .  (4.8) 

The  corresponding  real  part  X  'c(  ** ) =  H[  X  "c(  v  )  :  ]  is  then  obtained  as 

x'c(,,):IR  =  ‘  S  [  +  yf(  ** )  1  +  fg(-‘')  -  yf(-*'))j  •  (4.9) 

The  real  part  of  the  susceptibility  without  line  coupling  for  the  limiting  spectral  regions, 
given  by  Eqs.  (3. 8)  and  (3.13),  can  be  obtained  from  these  expressions  including  line 
coupling  by  setting  y  =  0. 

As  in  the  development  for  the  case  without  line  coupling,  we  examine  expressions 
for  two  limiting  spectral  regions,  Eqs.  (4.7)  and  (4.9),  and  postulate  the  following  form 
for  the  real  part  of  the  susceptibility  with  line  coupling 

X  'Cp(  "  )  =  S  tanh( 0  v  Q/2)  [  (-  g(v)  +  (  «  /  i»0)f(  »  )]  +  y[f(»>)  +  («/  ■'Q)g(i')] 

+  [-  g(-v)  +  (a  /  »'0)f(-  f)]  +  y[f(-  »)  +  ('  I  vj  g(-  v>]  j  .(4. 10a) 

Two  equivalent  forms  for  this  proposed  function  are  given  as 

X  ’cp(  v  )  -  S  tanh(  0  »Q/2)  [-(1  -  ya/v^giv)  +  (a  /  vQ  +  y)f(  v) 

-  (1  -  y«/i'0)g(-i')  +  (<*  /  •'g  +  y>f(- •  (4.10b) 

and 


x’cp(0  *  (8/  »  )  tanh(0  v  fl/2) 

l-  (•'-,,Q)  +  “2/i,Q  +  yo*’/»'0]  l(i»+K0)+o2/»'g-yo»/v0) 

[(  •>  -■'g)2  +«2]  [(  r  +  *'Q)2  +  «2] 


(4.10c) 
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The  agreement  of  the  proposed  function,  Eq.  (4.10),  with  expressions  for  the  lim¬ 
iting  spectral  regions  may  readily  be  established  by  allowing  0  v  q  «  1  for  the  millimeter 
region  and  by  taking  0  v  q  »  1,  a  /  0  and  y  -  1  for  the  infrared  region.  We  next 

examine  the  validity  of  this  real  part  by  taking  the  inverse  Hilbert  transform  and  com¬ 
paring  the  result  with  the  original  imaginary  part.  Using  Eq.  (4.10b),  we  derive 

X'cpO')  =H_1t  X  ’cp(  •')]  =  -  H[  x  ,cp(*’)] 

=  S  tanh(  0  v  g/2)  [[(1  -  y  a  /  vQ)f(  v)  +  (a  /  +  y)g(v )] 

-  [(1  -  y  “  /  *'0)f(-*')  +  (a  /  »Q  +  y)g(-  »)]  ]  .  (4.11) 

The  millimeter  limit  can  easily  be  checked  using  the  approximation,  tanh(  P  Vq/2)  =» 
PvqI2.  On  the  other  hand,  the  infrared  approximation,  tanh(  0  t>  q/2)  =»  1,  and 
a/t'Q  “0,  once  again  gives  the  correct  result.  Thus  the  appropriateness  of  the  pro¬ 
posed  real  part  given  by  Eq.  (4.10)  has  been  demonstrated. 

5.  LINE  WING  CONTRIBUTIONS 

The  question  arises  as  to  whether  the  real  part  of  the  susceptibility,  obtained 

through  the  Hilbert  transform,  depends  critically  on  the  line  wing  of  the  imaginary 

part.  This  issue  is  an  important  one  because  if  it  does,  then  the  use  of  the  impact 

line  shape  would  be  expected  to  introduce  a  significant  error  in  the  result  for  the  real 

part.  This  is  because  the  impact  line  shape  is  known  to  be  inappropriate  for  the  line 
22 

wings.  This  section  is  devoted  to  the  consideration  of  this  important  point. 

An  intuitive  approach  would  be  to  compare  a  line  shape  function  having  a  sub¬ 
stantially  different  form  for  the  line  wing  with  the  impact  result.  While  the  use  of  the 
Gaussian  function  for  X  "(  v)  would  serve  this  purpose,  the  lack  of  a  simple  closed 
form  expression  for  the  Hilbert  transform  of  the  Gaussian  suggests  the  desirability  of 
an  alternative  approach.  Let  us  consider  a  function  for  X”(v)  constructed  as  the 
difference  of  two  Lorentzian  functions  having  different  strengths  and  widths.  This 
approach  has  two  advantages:  the  Hilbert  transform  is  available  and  the  difference 
function  can  be  made  to  have  an  asymptotic  behavior  significantly  different  from  that 
of  a  single  Lorentz  function. 

22.  Clough,  S . A . ,  Kneizys,  F.X . ,  Davies ,  R .  ,  Gamache ,  R  . ,  and  Tipping,  R .  (1980) 
Theoretical  line  shape  for  H2O  vapor:  application  to  the  continuum,  in  Atmos¬ 
pheric  Water  Vapor,  A.  Deepak,  T.D.  Wilkerson,  and  L.H.  Ruhnke,  Eds., 
Academic  Press,  New  York,  NY,  pp  25-46. 


mwwww^wwiw 


For  clarity  of  the  presentation,  we  consider  the  simple  Lorentz  function  in  the 
infrared  region.  Let  L(  v )  and  L^(  v )  be  the  simple  Lorentz  and  the  modified  Lorentz 
profiles  centered  at  v  -  v  q  and  be  defined  by 

L(»)  =  f(«).  (5.1) 


Ld(v)  *  m1  f(»)  -  m2  f2(v),  (5.2) 

where  the  Lorentz  function  f(  v  )  has  a  line  width  <*  ,  and  f2(  v  )  in  Eq.  (5.2)  has  a  line 
width  a  The  integrated  strengths  of  the  two  profiles,  L(v)  and  L^O),  are  main¬ 
tained  through  the  relationship 


ral  '  m2  = 


Equation  (5.2)  may  be  rewritten  as 


(5.3) 


Ld<w)  ~T 


KitijO  -m2a2)(  V  -  !>0)  +  a  o  2(m1a  2-m2a  )] 


(5.4) 


[O'  ->'0)2+“2]  [O'  "  *,0) 2  +  <*  22] 


Setting 

mja  =  m2  “  2  ’  (5-5) 

the  v  dependence  of  the  numerator  of  Eq.  (5.4)  is  eliminated.  As  a  consequence,  the 

-  4  -  2 

new  line  function  asymptotically  approaches  zero  as  v  rather  than  v  of  the  original 
function.  Thus  a  simple  modification  has  provided  a  significant  difference  in  the 
asymptotic  behavior  of  the  line  shape. 

We  define  the  ratio  of  the  halfwidths  of  the  two  Lorentz  functions  in  Eq.  (5.2)  to 
be  r,  namely,  r  =  a  I  a  Then  Eq.  (5.5)  can  be  written  as 


Then  Eqs.  (5.3)  and  (5.6)  can  be  combined  to  give 


1  (1  -  r) 


(5.7) 
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r 


(5.8) 


J2  "  (1  -  r) 


Because  our  intention  is  to  alter  the  line  wing  without  appreciably  affecting  the  center, 


we  assume  the  ratio  r  to  be  small. 


Using  these  two  line  shape  functions  and  the  infrared  approximation  in  Eq.  (3.14), 
the  real  parts  x  '( v  )  and  x  'd(  v )  can  be  obtained  as 


x'(iO  =  -  S  g(v),  and 


x'  (v)  =  -  S  [ - g(»0 - g2(v)]. 

a  1  -  r  1  -  r 


(5.10) 


where  g(  v  )  and  v )  have  halfwidths  a  and  a 2  respectively.  After  some  manipula¬ 
tion,  Eq.  (5.10)  becomes 


x'H('')  »  X '( •*  )  I  1  +  r  ll  - 


(»  -•'0)2  +  r2  “22 


(.  -  ,0)2+  «22 


(5.11) 


Now  we  can  compare  this  with  the  original  impact  result.  Far  from  the  line  center,  for 
which  :  v  -  •'q:»  « 2 •  the  inn®r  bracket  vanishes  and  Eq.  (5.11)  reduces  to 


X'd(,,):far~  X*<  • 


(5.12) 


Thus,  the  important  result  is  obtained  that  the  real  part  of  the  susceptibility  far  from 
the  line  center  is  relatively  unaffected  by  the  change  in  the  asymptotic  property  of  the 
imaginary  part. 

On  the  other  hand,  near  the  line  center  where  :  v  -  v  g:«  a2’  leads 


xV^neBr-  (1  +  r>  x'^): 


(5.13) 


through  the  second  order  in  r.  This  small  change  in  magnitude  may  be  associated  with 
the  change  in  the  magnitude  near  the  line  center  of  the  imaginary  part,  namely,  for 


:  v  -  v  „:«  «2. 


xV»)!n*f”  <1+F)  X"(*0: 


near  ’ 


(5.14) 


'  ^  l-  .W/.VAvWaW. 


r»  »  »  •  »  •  » 


through  the  first  order  in  r.  Comparison  of  Eq.  (5.13)  with  Eq.  (5.14)  suggests  that 
the  slight  change  in  the  magnitude  of  the  real  part  near  the  line  center  is  attributable 
to  the  small  change  in  the  imaginary  part  near  the  line  center. 

To  provide  a  sense  of  the  magnitude  of  r,  let  us  consider  numerical  values  appro¬ 
priate  to  the  molecular  collision  broadening  problem.  We  take  a  to  be  0.1  cm  1  and  if 
we  assume  a  ^  to  be  10  cm  *,  we  obtain  a  value  for  r, 

0.1 

r  = - =  0.01  .  (5.15) 

10 

It  should  be  noted  that  for  small  values  of  r,  the  ratio  of  the  magnitudes  of  the  two 
real  functions,  x"^  (v0  ±  a2)/x"  ^v0  -  a2^  ’  ls  aPPr03dmately  one-half.  Thus,  when 
r  =  0.01,  the  modified  line  profile  is  about  one-half  the  value  of  the  Lorentz  profile 
10  cm'1  away  from  the  line  center.  The  magnitude  of  the  real  part  near  the  line 
center  is  increased  by  the  factor  1.01  (Eq.  (5.13)),  and  far  from  the  line  center, 

:v  -  v-:  ■»  10  cm  1 ,  the  real  part  is  unaffected  by  the  reduction  of  the  line 
wing  (Eq.  (5.12)). 

6.  SERIES  EXPRESSION  FOR  EXACT  REAL  PART 

We  now  address  the  evaluation  of  the  exact  real  part  of  the  susceptibility.  An 

exact  expression  for  X  '(  v )  in  an  infinite  series  form  is  derived  using  a  known  expan- 

23 

sion  of  the  tanh( . )  function  using  Lorentz  functions  is  given  by  Gradshteyn. 

A  1 

tanh(  »  x/2)  =  (4x/ »  )  2_ -  (6.1) 

k=l  x2  +  (2k  -  l)2 

We  substitute  this  in  Eq.  (4. 1)  to  obtain 


23.  Gradshteyn,  I.S.  and  Ryzhik,  l.M.  (1980)  Table  of  Integrals,  Series,  and 
Products  A.  Jeffrey  Trans.,  Academic  Press.  New  Yonc,  NY. 


where 


(2k  -  1)  ir 

e(k)  - - - - .  (6.3) 

This  result  provides  a  series  expansion  of  the  impact  line  shape  including  line  coup¬ 
ling.  From  this  expression  it  is  straight-forward  to  obtain  the  Hilbert  transform  of 
x"(v).  In  fact,  the  partial  fraction  expansion  of  the  summand  gives 


A  r  2A(k)* 

X"(«')=4S//3  2_  |_ -  +  [B(k)f(  v)  -  A(k)g(»/)1 


k=l  K2+e(k)2 


[B(k)f(  -v)  -  A(k)g(- v))]  ,  (6.4) 


with 


A(k) 


<»(  •'02+a2-e(k)2)  -  y  v  Q(,  v  a2+e(k) 2) 


D(k) 


(6.5a) 


B  (k)  = 


V  v02+“2+e(k)2)  +y°(<'02+“2-e(k)2) 


D(k) 


(6.5b) 


D(k)  =  (  *'02+or2-e(k)2)2  +  4  v  Q2e(k) 2  . 


(6.5c) 


Because  the  series  converges  absolutely,  the  Hilbert  transform  and  the  summation  can 
be  interchanged.  Thus  the  desired  result  is  obtained. 


X ’(  »)  =  H[  X"(v)] 


=  (4S/  0) 


E[ 

k=l 


2A(k)e(k) 
v2  +  e(k)2 


-  [B(k)g( v  )  +  A ( k ) f (  v )]  -  f B ( k)g(  -  v) 


-  fB(k)g(  -  v)  +  A(k)f(  -  v)] J 


(6.6) 
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This  series  expression  will  be  used  for  the  numerical  verification  of  the  accuracy  of  the 
proposed  real  part  X  '  (  v ) . 

For  ease  of  computation ,  we  rewrite  this  equation  in  the  form : 

"  2A(k)e(k) 

X'(v)  =  (4S//S)  V  - 

o  2 
k=1  lv  +  e(k)  ] 

oo 

-  (4S//J)  [f(»0  +  f(-iO]  £  A(k) 

k=l 

OO 

-  (AS  /  0  )  [gOO  +  g(-«>)l  £  B(k)  .  (6.7) 

k=l 

Although  the  first  summation  has  to  be  computed  for  each  wavenumber  value  v  ,  the 

second  and  third  summations  need  only  be  evaluated  once  for  a  given  transition  since 

they  are  independent  of  t> .  It  is  noted  that  the  convergence  rates  of  the  series  are 
3  2 

of  the  order  of  1/k  for  the  first  summation  and  of  1/k  for  the  other  two  summations. 
This  expression  provides  a  viable  procedure  for  calculating  X  ’( v  )  from  X  "(  v )  in 
the  impact  approximation  for  a  limited  number  of  transitions  at  a  limited  number  of 
wavenumber  values. 

The  exact  real  part  of  the  susceptibility  in  Eq.  (6.7)  and  the  proposed  function 
given  by  Eq.  (4.10)  have  been  compared  numerically  for  three  values  of  the  line  cen¬ 
ter,  v.  =  10,  410,  and  1500  cm  *,  which  represent  transitions  located  in  the  millimeter, 
U  _  1 
intermediate,  and  infrared  regions,  respectively.  The  choice  of  v  ^  =  410  cm  has 

been  made  because  it  provides  0  v  q/2  =*  1  for  room  temperature,  T  =  296  (K),  where 
neither  the  approximation  for  tanh(.)  of  Eq.  (3.5)  or  of  Eq.  (3.9)  is  appropriate. 
Typical  values  of  S  =  2.0  cm  1  and  a  =  0. 1  cm  *  have  been  used  for  computations 
throughout . 

Two  runs  have  been  performed  for  each  value  of  the  line  center:  once  for  the 

line  center  region,  v  -  -1  cm  1  to  v  n  +  1  cm  '  at  a  sampling  interval  of  0.01  cm  *, 

0  0  4-1 

and  the  second  for  wide  ranges  spanning  from  1  to  10  cm  for  the  millimeter  case  and 
from  10  to  10^  cm  1  for  other  cases  at  appropriate  sampling  intervals. 

The  line  coupling  effects  have  been  evaluated  by  considering  doublets  centered  at 
wavenumber  values  •'g±a>  providing  a  separation  2  a  =  0.2  cm  and  with  the  coupling 

coefficients  y.  =  ±0.1.  The  convergence  criterion  that  the  ratio  of  the  new  iteration 

*  -  fi 

term  to  the  accumulated  sum  be  less  than  a  prechosen  small  number  e  =  10  has  been 


utilized  for  terminating  the  first  summation  in  Eq.  (6.7),  and  the  numbers  of  iterations 
needed  for  convergence  have  been  recorded  for  each  run.  The  two  other  summations, 
evaluated  only  once  for  each  line,  have  been  computed  to  a  specified  maximum  number 
of  iterations  for  which  satisfactory  convergence  is  assured.  Additional  trial  runs  were 
also  performed  for  various  convergence  critera,  e  =  10  ,10  ,  and  10  ,  to  assess 

the  convergence  property  of  the  first  summation  of  Eq.  (6.7).  The  results  showed  20, 
40,  and  80  iterations  on  the  average  to  reach  the  convergence,  in  contrast  to  an  aver- 

_  g 

age  of  180  iterations  for  the  runs  with  e  =  10  .  This  indicates  that  the  series  expres¬ 

sion  obtained  for  the  exact  real  part  could  be  used  in  computational  application,  pro¬ 
vided  the  number  of  transitions  is  small. 

The  results  obtained  are  summarized  in  Figures  1  through  3  and  Table  1.  Figures 
1-A  through  1-C  represent  the  line  center  computations  for  millimeter,  intermediate, 
and  infrared  regions,  respectively,  for  a  single  line.  Each  figure  consists  of  three 
plots;  the  imaginary  part  X  "(  *),  the  proposed  real  part  X  '  (  v)  of  the  susceptibility, 
and  the  difference  between  the  proposed  and  the  exact  real  parts  x  '  <  v  )  -  x  '<  a ) . 
Figures  2-A  through  2-C  indicate  the  difference  over  wider  regions  for  the  single  line 
cases.  Figures  3-A  through  3-C  represent  computations  for  two  lines  with  line  coup¬ 
ling.  Table  1  summarizes  the  error  computations.  It  includes  the  number  of  terms 
required  for  the  convergence  of  the  first  summation  in  the  series  expression,  Eq.  (6.7), 
for  the  exact  real  Dart  of  the  susceptibility . 

7.  DISCUSSIONS  AND  CONCLUSIONS 

A  unified  and  analytic  approximate  real  part  of  the  susceptibility  over  the  entire 
millimeter  and  infrared  regions  that  is  consistent  with  the  Van  Vleck- Huber  impact  line 
shape  with  line  coupling  has  been  proposed.  Fundamental  to  our  derivation  is  the  as¬ 
sumption  of  superposition,  namely  that  the  complex  susceptibility  for  the  system  may  be 
obtained  from  the  contributions  to  the  susceptibility  from  the  individual  transitions. 

This  assumption  allows  us  to  apply  the  Hilbert  transform  and  its  inverse  to  the  real  and 
imaginary  parts  for  a  single  transition  to  establish  the  Kramers- Kronig  relationship  for 
the  system.  The  proposed  real  part  of  the  susceptibility  has  been  obtained  through 
careful  examination  of  two  limiting  cases  of  the  millimeter  and  infrared  spectral  regions 
for  which  the  real  parts  have  been  explicitly  computed  from  the  respective  limiting  forms 
of  the  imaginary  part.  The  inverse  Hilbert  transform  of  the  proposed  real  part  was  com¬ 
puted  to  demonstrate  consistency  with  the  original  function. 

Although  for  the  line  coupling  case,  the  zero-sum  rule  is  an  important  element  of 
the  spectral  property  of  both  the  real  and  imaginary  parts  of  the  susceptibility  for  the 
molecular  system,  its  utilization  has  not  been  necessary  in  the  development  of  the  pro¬ 
posed  function.  An  important  requirement  in  our  formulation,  necessary  to  maintain 


Figure  1A.  Imaginary  part  of  the  susceptibility  (X  "),  proposed  real  part  of  the 
susceptibility  (  x  ’p),  and  the  difference  between  the  proposed  and  actual  real 
part  (  X  ’p  -  X  ')  as  a  function  of  wavenumber  for  a  single  line  with  strength  of 
2  cm"1  at*10  cm"1  and  a  halfwidth  of  0.1  cm'1 
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Figure  IB.  Imaginary  part  of  the  susceptibility  (  X  ").  proposed  real  part  of  the 
susceptibility  (  X  ’p) ,  and  the  difference  between  the  proposed  and  actual  real 
part  ^X  'p  -  X  ')  as  a  function  of  wavenumber  for  a  single  line  with  strength  of 
2  cm' '  at  410  cm'  *  and  a  halfwidth  of  0.1  cm"  * 
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Figure  2A.  Imaginary  part  of  the  susceptibility  (  x  ").  magnitude  of  proposed  real 
part  of  the  susceptibility  (  X  'p),  and  the  difference  between  the  proposed  and 
actual  real  part  (  X  'p  -  X  ')  tfver  an  expanded  spectral  range  for  a  single  line  with 
strength  of  2  cm'1  at  10  cm'1  and  a  halfwidth  of  0. 1  cm"1 
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Figure  3B.  Imaginary  part  of  the  susceptibility  (  X  "),  proposed  real  part  of  the 
susceptibility  (  x  'p).  and  the  difference  between  the  proposed  and  actual  real 
part  (  x  'p  -  X  ')  as  a  function  of  wavenumber  for  two  coupled  lines  of  strength 
2  cm'*  centered  at  410  ±0.1  cm'l,  with  halfwidths  of  0.1  cm'*  and  coupling 
coefficients  of  ±0.1  respectively 
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Figure  3C.  Imaginary  part  of  the  susceptibility  (  x  "),  proposed  real  part  of  the 
susceptibility  (  X  'p),  and  the  difference  between  the  proposed  and  actual  real 
part  (  X  'p  *  X  ')  as  a  function  of  wavenumber  for  two  coupled  lines  of  strength 
2  cm’1  centered  at  1500  +0.1  cm'1,  with  halfwidths  of  0.1  cm'1  and  coupling 
coefficients  of +0.1  respectively 
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Table  1.  Summary  of  the  Accuracy  Test.  The  values  in  the  parentheses  are  modi¬ 
fied  percent  errors  where  a  single  extremum  percent  error  at  the  zero  crossing  is 
eliminated  from  each  case.  Other  values  are  Sj  -  2.0,  aj  -  0.1,  yj  =  -0.1,  y2  =  0.1, 
b'  =  1/412,  and  e  for  the  termination  criterion  is  10~6 


Vj/l/cm) 

Millimeter 

10.0 


9.9  &  10.1 

Intermediate 

410.0 

409.9  ft  410.1 

Infrared 

1500.0 


Range(l/cm) 

Average  Percent  Error 

No.  of  Iterat 

[9,  11] 

0.002 

50 

[0,  100] 

0.42 

50  -  51 

[0,  1500] 

22 

50  -  88 

[9,  11] 

0.004 

50 

[0,  100] 

0.46 

50  -  51 

[0,  1500] 

23 

50  -  88 

[409,  411] 

0.22(0.006) 

68 

[0,  1910] 

0.23(0.009) 

62  -  116 

[409,  411] 

0.15 

68 

[0,  1910] 

0. 11(0.010) 

62  -  116 

[1499,  1501] 

0.97(0.006) 

177 

[0,  3000] 

0.65(0.003) 

136  -  213 

[1499,  1501] 

0.47(0.035) 

177 

[0,  3000] 

0.29(0.003) 

177  -  213 

consistency  of  the  real  and  imaginary  parts,  has  been  the  existence  of  a  complete  pro¬ 
file  containing  both  positive  and  negative  oscillators.  Constant  biases  in  these  oscil¬ 
lators  cancel  each  other  and,  therefore,  the  resulting  function  vanishes  at  infinity. 

The  Hilbert  transform  of  the  resulting  function  is  thus  well-defined. 

The  effect  of  the  line  wing  of  the  imaginary  part  on  the  Hilbert  transform  has  been 

assessed  using  a  profile  composed  of  the  difference  of  two  Lorentz  functions,  for  which 

-  4  -  2 

the  line  wing  decays  as  v  instead  of  v  for  the  simple  Lorentz  shape.  The  real  part 
in  the  region  far  from  line  center,  obtained  from  the  modified  profile,  has  been  shown 


.’VA".'  • 


CTO 


to  have  the  same  asymptotic  behavior  as  that  of  the  simple  Lorentz  case.  This  demon¬ 
strates  the  decreased  sensitivity  of  the  Hilbert  transform  to  the  line  wing  of  the  imagin¬ 
ary  part  and  supports  the  use  of  the  impact  line  shape. 

For  evaluation  of  the  validity  of  our  approximate  function,  the  Hilbert  transform  of 
the  imaginary  part  of  the  susceptibility  has  been  obtained  in  the  form  of  an  infinite 
series  by  exploiting  the  Lorentz  expansion  of  the  tanh(.)  function  and  the  known  Hilbert 
transform  of  the  Lorentz  function.  The  resulting  series  has  been  successfully  used  in 
numerical  verification  of  the  proposed  approximate  real  part.  Excellent  agreement  has 
been  demonstrated  in  the  vicinity  of  the  line  center  and  generally  good  agreement  over 
wider  ranges.  However,  the  gradual  deterioration  in  the  percentage  accuracy  far  from 
line  center  in  the  millimeter  case  may  be  a  significant  problem  for  certain  applications. 

Several  approaches  could  be  considered  to  improve  the  validity  of  the  proposed 
function.  A  fundamental  reason  for  the  error  far  from  line  center  associated  with  the 
millimeter  case  is  attributable  to  the  inadequate  approximation  of  the  tanhf.)  function 
over  a  broad  spectral  range.  Consequently,  for  a  broad  spectral  range  it  is  expected 
that  an  improved  approximation  for  the  tanh(.)  function  would  be  required  in  order  to 
obtain  improved  results.  Alternatively,  empirical  corrections  could  be  developed  based 
on  the  functional  characteristics  of  the  error.  Of  particular  significance  would  be  the 
introduction  of  a  small  empirical  shift  in  the  line  center  to  substantially  reduce  the  per 
centage  error  in  the  line  center  region  at  which  the  real  part  of  the  susceptibility  goes 
through  zero. 

Finally,  it  is  noted  that  the  proposed  function  can  be  used  in  conjunction  with  the 

line  parameters  over  the  entire  frequency  range.  Thus  the  vast  data  base  such  as  the 

line  parameter  compilation  can  be  used  for  the  evaluation  of  the  real  part  of  the 

susceptibility.  The  nrnpnsoil  function  can  be  readily  incorporated  into  line-by- 

24*26 

line  computational  code. 


24.  Clough,  S.A.,  Kneizys,  F.X.,  Rothman.  L.S.,  and  Gallery,  W.O.  (1981)  Atmos 

pheric  spectral  transmittance  and  radiance:  FASCOD1B.  SP1E  VoL27J  Atmos¬ 
pheric  Transmission,  pp  152-166.  ~  ~ 

25.  Rothman,  L.S.,  et  al.  (1983)  AFGL  trace  gas  compilation:  1982  version,  Appl. 

Opt.  ,  22,  No.  11:1616-1627.  - 

26.  Rothman,  L.S.,  et  al.  (  1983)  AFGL  atmospheric  absorption  line  parameters  compi¬ 
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In  this  appendix,  a  summary  of  the  classical  relationships  among  the  susceptibil¬ 
ity,  dielectric  constant,  and  index  of  refraction  are  presented.  The  definitions  of  the 

A1 

parameters  and  the  derivation  of  the  wave  equation  follow  Abraham  and  Becker. 

For  a  homogeneous,  isotropic  dielectric  in  which  the  current  density  and  charge 
density  have  been  assumed  to  be  zero.  Maxwell's  equations  are 

1  3D 


v  x  H  =  — 
c 

3 1 

(Al) 

1 

vx  E  = - 

c 

d§ 

at 

(A2) 

v-  5  =  o. 

(A3) 

v-  B  =  0, 

(A4) 

where 

H  =  magnetic  field  strength 

E  =  electric  field  strength 

D  =  dielectric  displacement 
B  =  magnetic  induction 

and  c  is  the  velocity  of  light . 

We  also  have  the  further  definitions,  for  non -ferromagnetic  substances,  following 


A1 

Abraham  and  Becker  that 

D  =  E  +  4  »P,  (A5) 

where  P  is  the  dielectric  polarization  and  that 

D  =  e  E,  (A6) 

B  =  m  H,  (A7) 

P=  XE,  (A8) 


Al.  Abraham,  M.  and  Becker,  R,  (1949)  The  Classical  Theory  of  Electricity  and  Mag¬ 
netism.  Hafner  Publishing  Company,  New  York,  NY. 


e  =  dielectric  constant 
=  permeability 
X  -  dielectric  susceptibility 


Substituting  Eqs.  (A6)  and  (A8)  into  Eq.  (A5),  and  neglecting  any  internal  electric 
fields  gives 


*  =  1  +  4  *  x  .  (A9) 

A 1 

As  shown  in  Abraham  and  Becker,  the  wave  equation,  derived  from  Maxwell's 
equations  and  Eqs.  (A6)  and  (A7),  is  given  by 


2  2 
c  at 

with  a  similar  equation  for  the  magnetic  field  strength  H. 

The  velocity  of  propagation  of  the  electromagnetic  wave,  v,  is  defined  as 

v  =  c  IW*  , 

and  the  refractive  index  of  the  medium ,  n ,  as 
n2  =  (c/v)2  =  Me  . 

For  materials  with,  m  *»  1,  Eq.  (A  12)  becomes 


(A  10) 


(All) 


(A12) 


n2  =  «  ,  (A13) 

and  we  have,  by  combining  Eqs.  (A9)  and  (A13) 

4  ir  X  =  n2  -  1,  (A14) 

with  n,  e  ,  and  X  real. 

For  a  plane  wave  of  amplitude  Eg  and  wavenumber  value  v  =  w/2  it  c,  propagating 
in  the  z-direction,  the  solution  to  Eq.  (A  10)  is 

39 


I 


Ex  =  Eg  exp[i2  ir  v  (ct  -  nz)] ,  (A15) 

with  Ez  =  Ey  =  0,  and  a  similar  equation  for  the  magnetic  field  Hy  of  amplitude 

Hy  =  Hg  exp{i2  tv  (ct  -  nz)] ,  (A16) 

with  Hx  =  Hz  =  0.  Substituting  Eqs.  (A15)  and  (A16)  into  Eq.  (A2)  gives  the  relation¬ 
ship  between  the  field  amplitudes. 


H0=TV  (A1' 

The  energy  transport  per  unit  area  per  unit  time  is  given  by  the  Poynting  vector  N, 
defined  by 


N  = -  (E  x  H).  (A  18) 

4  t 

and  for  the  plane  wave  solution,  the  cycle  averaged  component  of  the  Poynting  vector 


1  c  n  , 

<  N  >=-(—)  (— )  E  . 
z  2  At  P  u 

For  the  general  case  where  the  dielectric  constant  is  complex,  that  is 


«  =  e  '  -  i  e  ” , 


the  index  of  refraction  and  the  susceptibility  also  become  complex  and  are  defined  as 


n  =  n'  -  in", 


X  =  X  1  -  i  X  ", 


Equation  (A  12)  is  still  valid  with  the  complex  quantities,  and  the  plane  wave  solution  is 


Ex  =  Eg  exp[i2  t  v  (ct  -  nz)] 


=  Eg  exp]-2irn"  vz]  exp[i2  *  v  (ct  -  n'z)] , 


( A23) 


for  the  electric  vector,  and 


Hy  =  H0exp[-2»n"  rz]exp[i2  i  v  (ct  -  n'z)l , 


(A24) 


for  the  magnetic  vector,  with 


n'-in" 


Ho=— Eo=  — 


E0- 


CA25) 


An  alternate  form  of  Eq.  (A25)  explicitly  in  terms  of  the  phase  angle  between  the 
electric  and  magnetic  amplitudes  is  given  by 


|n|  -i# 

H0  =“e  E0’ 


(A26) 


where  tan  0  =  n"/n'  and  |n|  =  [n'2+n"2] 

If  n  *  1,  Eq.  (A  13)  can  be  written  explicitly  in  terms  of  its  real  and  imaginary 

parts 


or 


From  Eqs,  (A9)  and  (A29),  for  the  susceptibility  we  obtain 


4  *  X  '  =  n’2  -  1-n"2, 

(A  30) 

4*  X  "  =  2n'n". 

(A31) 

For  a  complex  index  of  refraction  and  dielectric  constant,  Eqs.  (A23)  and  (A24) 
are  damped  plane  wave  solutions  with  the  real  part  of  the  cycle  averaged  Poynting  vec¬ 
tor  given  by 
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K 

(n*  -  i  n")  =  c  T  -  i  €  ", 

(A27) 

5$ 

> r  J\ 

i 

?  2 
c  '  =  ir  -  n”  , 

(A28) 

4 

e  "  =  2  n’n". 

(A29) 
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<nz> 


1  c  n'  , 

—  ( - )  C— )  E^exp  [-4*n"*z], 

2  4  it  M 


(A  32) 


which  defines  an  absorption  coefficient ,  «  ,  as 


k  =  4in"v. 
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Appendix  B 


Hilbert  Transform  of  Lorentz  Type  Functions 
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We  consider  the  Hilbert  transform  of  the  following  two  Lorentz  type  functions  that 

appear  repeatedly  in  the  main  body  of  this  report.  We  note  that  this  appendix  follows 
B 1 

Brown  closely. 


f  ( v  )  = 


(o  /  n  ) 

[(<"  »0)2  +  “21 


g(  v  )  =  - 


[(  p  -  pQ)2  +  a2] 


The  Hilbert  transform  of  functions  f(-  v)  and  g(-  v)  which  also  appear  in  the  text 
can  be  obtained  from  those  of  the  above  functions  as  shown  later. 

Consider  the  Hilbert  transform  F(»')  of  (Bl)  first. 


OO 

/(<*/») 
- - 

[(U  -  v  )  (u  -  vnr  +  a 


Using  the  change  of  variable 


we  can  rewrite  the  integral  in  Eq.  (B3)  as 


F(«0  =  (1/  *  a) 


[ (t  -  x)  (t‘  +  1)] 


where  x  =  (*»-»>  q  )/<*.  The  Cauchy  principal  value  of  the  integral  in  Eq.  (B5)  can 
be  obtained  by  expanding  the  integrand  using  the  partial  fraction  as 


Bl.  Brown,  D.R.  (1984)  Dispersion  profiles  for  near  millimeter  wave  oxygen  refractiv- 
ity,  Private  communication. 


[  *  a  (x  +1)]  (t  +1)  (t  +1) 


=  -gO)- 


1)] 


On  the  other  hand,  the  change  of  the  dummy  variable  in  the  definition  of  the  Hilbert 
transform  from  u  to  -u  gives  us  the  following  property  for  a  general  function  0  . 


H[0(->O]  =  -  H[#0)]: 


Therefore  the  Hilbert  transform  of  f(-  v)  is 


H[f(-  r)J  =  -  F  (-!>>  =  g(-*). 


For  H[g(  k)]  ,  we  may  either  use  the  symmetry  of  the  Hilbert  transform0  and  the 
above  result  or  follow  a  similar  derivation.  Here  we  will  go  through  the  latter  path. 
After  the  same  change  of  variable,  we  arrive  at 


GOO  =  H[gO)l 


=  (1/  *  «) 


[(t  -  x)  (t4  +  1)] 


where  x  is  defined  above.  Following  the  F(  «0  case,  we  have 


GO)  = 


[  jr  2a 


1  r  x  xt  1 

(x2+l)]  J  (tx)  (t2+l)  (t2+l 


=  f  ( » ) . 


t»«(x‘+  1)] 


B2.  Bracewell,  R.N.  (1978)  The  Fourier  Transform  and  Its  Applications,  Second  Edi¬ 
tion,  McGraw-Hill,  New  York,  NY. 
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